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The recent Letter [l[ describes full quantum simula- 
tions of a dark soliton in a Bose-Einstein condensate in a 
regime where the system cannot be described by the per- 
turbative approach. The authors argue, based on the fill- 
ing in of the two-point correlator g'-^^ that a photograph 
of a condensate would reveal a smooth atomic density 
without any localized dark soliton. This is in contrast to 
the perturbative regime where a photograph would show 
a dark soliton with a random position While we ad- 
mire the quantum simulations and other results in [ij , we 
think that their conclusion about the outcome of a single 
experiment is not justified by this property of .g'^^ 

As the question of whether and when dark solitons fill 
in in single realisations has been a recurring and some- 
times confusing one in the field, it is important to pin 
down what conclusions can or cannot be made. With 
this aim, we provide the following counterexample in the 
non-perturbative regime (see also Q). 

Let 4>q(x) cx tanh[(a; — q)/^] be a standard conden- 
sate wave function with a dark soliton at q. With 


dg = j dx (l)*{x)'^{x), a state (aj)^ |0) is a condensate 
with a soliton at q, and let the iV-particle state be a 
superposition of condensates with different q 
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where 4'o{q) defines the superposition. 

After measurements of n atomic positions xi, a;„ the 
state ([T]) collapses to a conditional state 
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where tpn{q) = 4>q{xn)--<j)q{xi)'4jQ{q)- The (n + l)-st mea- 
surement will find a particle at Xn+i with a probability 

Pn+l{Xn+l) OC (V'n|^'^(a;n+l)*(a;ri+l)|V'n)- 

This is equivalent to simultaneous measurement of all Xi . 

Using the methods of we simulated measurement 
of all N = 5000 particles on a lattice of 31 sites, where 
x,q G { — 15,15}, assuming a soliton width ^ — 1.5, and 
a delocalized (uniform) superposition tpo{q) cc 1 that is 
non-perturbatively wider than the soliton width. The in- 
set in the figure shows the ensemble average particle den- 
sity pi (x) and the main figure shows a generic histogram 
of particle positions xi, ..jXn measured in a single real- 
ization. Each single realization of the experiment finds a 
soliton localized at some definite but random q. 
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FIG. 1: (Color online) Histogram of measured atom positions 
in a single experiment from ([l]). Inset: single particle density 
Pi{x) (dashed black line) and g2(x) (solid red line), normal- 
ized so that J g2{x)dx = 1. 


What about the two-point correlator g2{x) — 
(■0o|^'''(O)\^'l'(a;)^(a;)^(O)|^Ao) that is analyzed in [H? 
This turns out nearly uniform (see the inset). Hence, 
the "filled in (72" ^ "filled in soliton" line of reasoning 
is clearly incorrect. Of course, the soliton simulated in 
[l| may still be greying, but the point here is that one 
cannot answer such a question by analyzing 52(2;). 

Our example demonstrates that, in some cases, a low 
order correlator like g2{x) is insufficient to draw conclu- 
sions on the outcome of a single experiment, and that 
the soliton in a BEG is such a case. Here, 52 is equal 
to our P2{x2) after the first particle was measured at 
xi = 0. A measurement of only one particle is not 
enough to collapse the soliton position. If we want to 
infer the soliton position from a histogram of particle po- 
sitions, then the number of measured particles must be 
large enough to provide a histogram with a well-resolved 
soliton notch. This is a fundamental requirement for an 
accurate Bayesian inference of the soliton position. 
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